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ABSTRACT 

This paper proposes and discusses a test of the chemical composition of ultra-high energy cosmic 
rays that relies on the anisotropy patterns measured as a function of energy. In particular, we show 
that if one records an anisotropy signal produced by heavy nuclei of charge Z above an energy I?thr, 
one should record an even stronger (possibly much stronger) anisotropy at energies > Ethi/Z due to 
the proton component that is expected to be associated with the sources of the heavy nuclei. This 
conclusion remains robust with respect to the parameters characterizing the sources and it does not 
depend at all on the modelling of astrophysical magnetic fields. As a concrete example, we apply this 
test to the most recent data of the Pierre Auger Observatory. Assuming that the anisotropy reported 
above 55 EeV is not a statistical accident, and that no significant anisotropy has been observed at 
energies < 10 EeV, we show that the apparent clustering toward Cen A cannot be attributed to heavy 
nuclei. Similar conclusions are drawn regarding the apparent excess correlation with nearby active 
galactic nuclei. We then discuss a robust lower bound to the magnetic luminosity that a source must 
possess in order to be able to accelerate particles of charge Z up to 100 EeV, Lb > 10 45 Z~ 2 erg/s. 
Using this bound in conjunction with the above conclusions, we argue that the current PAO data 
does not support the model of cosmic ray origin in active radio-quiet or even radio-loud galaxies. 
Finally, we demonstrate that the apparent clustering in the direction of Cen A can be explained by 
the contribution of the last few gamma-ray bursts or magnetars in the host galaxy thanks to the 
scattering of the cosmic rays on the magnetized lobes. 
Subject headings: cosmic rays 
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1. INTRODUCTION 

The sources of ultra-high energy cosmic rays have re- 
mained elusive in spite of the enormous progress reached 
on the experimental side, with present day detectors 
reaching apertures > 10 000 km 2 sr yr. The differen- 
tial energy spectrum, the chemical composition and the 
distribution of arrival directions on the sky are as many 
clues to the nature of the source. There is now a consen- 
sus on the existence of the GZK cut-off (Greisen 1966; 
Zatsepin & Kuz'min 1966), which has been observed by 
two different experiments (Abbasi et al. 2008, Abraham 
et al. 2008a). The existence of this GZK cut-off puts 
on solid ground the models which attribute the origin of 
ultra-high energy cosmic rays to powerful astrophysical 
objects distributed on cosmological scales, such as pow- 
erful radio-galaxies (Rachen & Biermann 1993), gamma- 
ray bursts (Milgrom & Usov 1995, Vietri 1995, Waxman 
1995), or magnetars (Arons 2003). 

Experimental results on the chemical composition and 
anisotropies are however far more confusing. While the 
most recent analysis of HiRes data points to a proton 
composition above the ankle (Abbasi et al. 2005, Belz 
2009), the fluorescence data of the Pierre Auger Observa- 
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tory (PAO) rather indicates an increasingly heavier com- 
position above this energy, with the last data points at 
~ 30 — 50 EeV close to expectations for iron (Unger et 
al. 2007, Wahlberg et al. 2009). Regarding the distri- 
bution of arrival directions on the sky, there exist vari- 
ous contradictory claims, see for instance the correlation 
with the super-galactic plane reported by Stanev et al. 
(1995), contradicted by the results of the AGASA ex- 
periment (Takeda et al. 1999); the claim for a correla- 
tion with BL Lac objects in Tinyakov & Tkachev (2001, 
2002), Gorbunov et al. (2002, 2004), which has been 
debated (Evans, Ferrer & Sarkar 2003, 2004; Tinyakov 
& Tkachev 2004); as well as the possible detections of 
multiplcts in various datasets, e.g. Takeda et al. (1999), 
Uchihori et al. (2000), Farrar, Berlind k Hogg (2006), 
the significance of which is questioned in Abbasi et al. 
(2004) and Finley & Westerhoff (2004); or, finally, the 
recent data of the PAO, which reveal a correlation with 
nearby active galactic nuclei (AGN), the maximal signal 
being obtained for a search radius of 3.1° around AGN 
located closer than 75 Mpc and for energies above 57 EeV 
(Abraham et al. 2007, 2008b). 

In general, these issues of anisotropies and chemical 
composition are discussed separately. However, as we 
argue in the present paper, one can construct a powerful 
test of the chemical composition by using the anisotropy 
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signal at various energies. This becomes particularly in- 
teresting when one notes that measurements of chemi- 
cal composition are relatively sensitive to the details of 
the shower reconstruction, in particular to the extrapo- 
lation of hadronic models at energies beyond those cur- 
rently tested in accelerators (see Wibig 2008, Ulrich et 
al. 2009a, b and Wibig 2009 for recent discussions of this 
issue). Moreover, measurements of chemical composition 
rely on the use of fluorescence detectors whose duty cycle 
is rather low and cannot be made event by event, so that 
they are limited by statistics. 

The essence of the test proposed in this paper is to 
exploit the fact that a source in the sky emitting heavy 
nuclei of charge Z at an energy E is expected to produce 
a similar anisotropy pattern at energies EjZ via the pro- 
ton component which is expected to be associated with 
the same source. Our expectation for the existence of a 
proton component relies on the theoretical expectation, 
that if protons are present in the plasma in which par- 
ticle acceleration takes place, then they should be accel- 
erated along side with the heavy nuclei. As discussed in 
some detail at the end of § a source accelerating heavy 
nuclei to energy E is quite generally expected to accel- 
erate protons to energy EjZ. Although the isotropic 
background "noise" increases at lower energies, we show 
in § [2] that the signal-to-noise ratio for the detection 
of the anisotropy also increases. This prediction does 
not depend at all on the modelling of astrophysical mag- 
netic fields as it only relies on the property that protons 
of energy EjZ follow the same path in the intervening 
magnetic fields and produce the same angular image as 
heavy nuclei of charge Z and energy E. The proposed 
test is discussed in more detail in Section O When ap- 
plied to the most recent data of the PAO (Section [3|), it 
shows that the signal that is responsible for the appar- 
ent anisotropy pattern at energies > 55 EeV must not be 
heavy but light nuclei, provided this anisotropy is con- 
firmed by future data, and provided the PAO does not 
see evidence for anisotropy at lower energies, as seems to 
be the case (Abraham et al. 2007). 

We discuss the implications of these results in the last 
two Sections. In Section [4] we argue that local AGN (in- 
cluding FR I radio-galaxies) do not possess the power re- 
quired to accelerate protons to ultra-high energies. When 
combined with the previous conclusions, this leads to the 
conclusion that the current PAO data do not support 
AGN as sources of the highest energy cosmic rays. Fi- 
nally, Section [5] summarizes the findings and concludes 
that the current data on composition and anisotropy 
suggest one of the following: (i) the shower modelling 
or, what would be more interesting, the hadronic the- 
oretical models of shower reconstruction are in error at 
high energy; (ii) the composition switches abruptly from 
heavy to light above 30 — 50 EeV; (iii) the source injects 
primarily heavy nuclei (which seems unlikely); (iv) the 
anisotropy seen by the PAO is a statistical artefact. 

2. TESTING THE CHEMICAL COMPOSITION OF 
ULTRA-HIGH ENERGY COSMIC RAYS WITH 
ANISOTROPY DATA 

As argued below, one can use the results of searches for 
anisotropy at various energies in order to constrain the 
chemical composition of ultra-high energy cosmic rays. 
The basic claim is the following: if one detects anisotropy 



above some energy E t hr, but not below, then the chem- 
ical composition is most likely light above E t hr, because 
one would have otherwise observed a similar anisotropy 
pattern at an energy Ethi/Z, with Z the assumed average 
charge of the cosmic rays at energies > Ethr . 

Consider a region of the sky, of angular size 9 and solid 
angle AVt = 2w(l — cos 8), showing an excess of particles 
above isotropic expectations above some energy thresh- 
old E t \ lr . The most natural interpretation is to infer the 
existence of one or more sources in this direction, the 
image of which is spread over Afi, because of smear- 
ing by 9 in the intervening magnetic fields, or because 
Afi subtends the source distribution. For the sake of 
the argument, we assume that the spectrum produced 
by the source consists of protons and of heavy nuclei 
of charge Z. In all known source models, one expects 
that the ratio of the elemental spectra q p (E)/qz(E) > 1 
and more generally q p (E) / 'qz(E) 3> 1 at a given energy 
E -c E max (p), where E max (p) denotes the maximal en- 
ergy at the source for protons. For instance, the com- 
position ratio of protons to iron peak elements (Z > 17) 
in the Galactic cosmic ray spectrum, which is roughly 
consistent with the solar abundance when compared at a 
given energy per nucleon, implies a source composition at 
a given energy q p : qz>n — 1 : 0.06 (as taken from the 
recent ATIC-2 data, Panov et al. 2006). The only energy 
regime in which one may obtain q p (E)/qz{E) < 1 is 
E max (p) < E < Bmax(^), with E max (Z) the maximal 
energy for the nuclei of charge Z (see also further below) . 
We also assume here that the maximal energy scales as 
Z; as argued further below, this is a conservative choice. 
Now, let us ask what would be seen if the composition 
above E t h r were heavy. Consider the quantity 

y (~>e )- AN{> Et ^ 

2^Z\> Ethr) — 



v/Mso(> E thl ) ' (1) 
which characterizes the signal-to-noise ratio of the 
anisotropy signal: AN(> Ethr) represents the excess 
number of events over isotropic expectations, i.e. the 
difference between the total number of events observed 
above -Ethr within Afl and N- lso (> Ethr), the number 
of events expected from the isotropic background. In 
the above model, AN(> Ethr) corresponds on average 
to the number of events produced by the sources in- 
side Afi. A detection of anisotropy thus corresponds 

tO S Z (> Ethr) > 1. 

One can now compute a similar quantity, but corre- 
sponding to the proton component at energies > E t hi/Z, 
denoted S p (> E t[i r/Z). These protons have a same 
magnetic rigidity than the heavy component above E t h r , 
therefore they follow the same path in the Galactic and 
extra-galactic magnetic fields, and consequently they are 
to produce a similar anisotropy pattern in the sky, up 
to the composition ratio of the elemental spectra and 
up to the level of background noise associated with 
the isotropic component. Assuming that the elemen- 
tal spectra at injection have the same power-law form 
<7i(E) oc E~ s at E -c E max (Z l ), and that the total 
number of events observed scales as E 1 ~ Sobs , with s bs 
the slope of the observed all-sky differential spectrum, 
one finds 

q p (E t hr/Z) 



S p (> E thT /Z) = S z (> E, 



thr 



qz(Ethi/Z) 



_£S-(sobs + l)/2 

(2) 
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The fudge factor ai oss > 1 incorporates the effect of 
energy losses; it is discussed in more detail in the follow- 
ing. At energies above the ankle, one has s b s — 2.7 
(Abbasi et al. 2008, Abraham et al. 2008a). There- 
fore, if the source emits a power-law with s = 2, one 
finds that the signal-to-noise ratio for detecting the 
anisotropy is larger at Ethr/Z than at E t hr by a factor 
aioss Z°- 2 q p (E thl ./Z)/q z (E thr /Z), generally expected 
to be significantly larger than unity! Clearly, the softer 
the source spectrum, the larger the gain in signal-to-noise 
ratio with decreasing Z. 

The factor ai oss accounts for the difference between the 
propagated and injected spectra. Given that the energy 
loss distance for protons of energy E t hr/Z is much larger 
than the energy loss distance for heavy nuclei of energy 
E t hr, if Ethr ^ 10 19 eV, one expects ai oss > 1- Indeed, 
regarding the (primary) proton component, energy losses 
can be safely neglected at E tri i/Z energies, since the at- 
tenuation length is of order of a Gpc or more at EeV en- 
ergies, which is much larger than the depth up to which 
which anisotropies can be produced, this latter being of 
order of the inhomogeneity scale of the large scale struc- 
ture ~ 100 Mpc or of the distance to the closest source, 
whichever is larger. Therefore, ai oss approximately cor- 
responds to the reciprocal of the attenuation factor for 
heavy nuclei of charge Z and energy > E t h r propagat- 
ing over the distance scale I to the sources inside Af2. 
The attenuation factor is understood as the ratio of the 
flux above > -Ethr after propagation on a distance I to 
the injected flux above this energy. For the case of iron 
primaries, this attenuation factor is 0.7 above 55 EeV 
at a distance of 10 Mpc, becoming 0.4 at a distance of 
100 Mpc (counting in the propagated flux all secondary 
nuclei with Z > 17). For oxygen, the attenuation factor 
is 0.5 at 10 Mpc, becoming 0.1 at 100 Mpc (counting 
in all C, N and O secondary nuclei); see Bertone et al. 
(2002) for an illustration of the relative fragility of ultra- 
high energy CNO nuclei vs iron peak elements. 

Actually, this enhancement factor ai oss should be fur- 
ther increased by the number of secondary protons with 
energy > Ethr/Z produced by the photo-disintegration 
of the primary heavy nuclei with energy > Ethr- This 
number is not easy to quantify but it can be significant. 
For instance, recent simulations by Aloisio, Berezinsky & 
Gazizov (2008) indicate that the signal should be domi- 
nated by protons at energies > 1 — 2 x 10 19 eV even if the 
source composition were strongly enriched in iron nuclei 
[assuming of course that E max (p) > 3 x 10 19 eV]. Al- 
lard et al. (2008) also conclude that heavy nuclei could 
be found in abundance at the detector only if the com- 
position were essentially dominated by Fe group nuclei 
at energies beyond 10 20 eV. As discussed in this latter 
reference, and above, the generic case that gives rise to 
a heavy composition at ultra-high energies is when the 
proton maximal energy at the source is smaller than the 
GZK cut-off. In the following, we adopt a conservative 
point of view and neglect the effect of energy losses, i.e. 
we set aioss = 1- 

Let us now discuss the scaling of E max with charge. 
The maximal energy is given by the comparison of the 
acceleration timescale (which depends on rigidity E/Z) 
with the minimum of the escape timescale (which also 
depends on E/Z), the age of the source or dynamical 
timescale (which do not depend on E, Z) and the en- 



ergy losses timescale (which depends on E, Z and mass 
number ^4). If the latter timescale does not provide the 
dominant limitation to the acceleration mechanism, then 
the maximal energy scales as the rigidity, as assumed 
above. If energy losses are dominant however, one ex- 
pects E max (Z) < ZE max (p), as explained in the follow- 
ing. Then the above argument would remain valid, as 
it only requires that protons exist at energies > E t hr/Z, 
i.e. that E max (p) > Ethr/Z- Note that in the more 
extreme case in which protons and heavy nuclei had sim- 
ilar maximal energies, the composition at ultra-high en- 
ergies > Ethr would be dominated by protons, and there- 
fore the anisotropy would be due to those protons since 
they have a higher rigidity than the heavy nuclei of a 
same energy. For photo-hadronic interactions, the en- 
ergy loss timescales for protons and iron-peak nuclei are 
known to be comparable, while that of intermediate el- 
ements is generally smaller; in that case, one would ex- 
pect E max (p) ~ E max (Fe) > E max (Z ~ 10). Regard- 
ing synchrotron losses, the energy loss timescale scales as 
(Z/A^E- 1 , therefore E max cx Z l / 2 {Z/A)- 2 if the ac- 
celeration timescale i acc cx E/Z and synchrotron losses 
dominate at the highest energies. In that case, the scal- 
ing E max (Z)/E max (p) cx Z remains correct to within a 
factor of order unity. The above conclusions agree with 
the recent simulations of Allard & Protheroe (2009). 

Finally, accounting for a more complex chemical com- 
position with more than two species would not modify 
significantly the results shown in Eq. (TT]) as long as the 
anisotropy above E tri r is indeed due to the species of 
charge Z . In practice, it will be useful to test the compo- 
sition above E t h r by checking the anisotropies at various 
Ethr/Zi, with Zi — 2, 8, 14, 26 representative of the most 
abundant elements. 

In the following, we provide specific examples of the 
above test applied to the most recent data of the PAO. 
We will show in particular that the above argument re- 
mains valid when one considers the all-sky anisotropy 
pattern due to an ensemble of close-by sources. 

3. APPLICATION TO THE PIERRE AUGER DATA 

Abraham et al. (2007, 2008b) have reported a corre- 
lation of 20/27 arrival directions of the highest energy 
events (E > 5.7 x 10 19 eV) with nearby (d < 75 Mpc) 
active galactic nuclei. However, the typical model of 
ultra-high energy cosmic ray origin in AGN refers to 
strongly beamed Fanaroff-Riley II (FR II) sources with 
giant radio lobes (e.g. Rachen & Biermann 1993), while 
19 out of these 20 correlating AGN in the Pierre Auger 
dataset belong to the Seyfert or LINER class, only one 
being a Fanaroff-Riley I (FR I) radio-galaxy (Cen A). 
These correlating AGNs are thus more likely to be trac- 
ers of the local large scale structure in which the actual 
sources of ultra-high energy cosmic rays camouflage and 
indeed, the arrival directions are compatible with a dis- 
tribution of sources following the large scale structure 
(Kashti k, Waxman 2008, Ghisellini et al. 2008, Zaw, 
Farrar & Greene 2009). The extended dataset up to 
2009 of the PAO reveals a significantly weaker correla- 
tion, with only 26 out of 58 events above 55 EeV lying 
within 3.1° of an AGN located within 75 Mpc (Hague et 
al. 2009). This level of correlation is compatible with 
that expected for sources distributed as the large scale 
structure (see discussion in Kotera & Lemoine 2008 and 
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Aublin et al. 2009 for a detailed discussion of the cor- 
relation with cosmological catalogs). It is important to 
note that the departure from isotropy lies at the 1% level. 
Therefore one should be careful not to overinterpret this 
excess. In the following, we use these PAO results as a 
concrete example for the test developed in the previous 
section. 

Interestingly, the PAO data up to 2007 has also shown 
evidence for clustering of events in the region around 
Cen A, with as many as 9 events within 20° (Gorbunov 
et al. 2008), as well as correlation of 8 out of 27 events 
with a sample of radio-galaxies (Moskalenko et al. 2008; 
Nagar & Matulich 2008). This has triggered a surge of 
interest in models of ultra-high energy cosmic ray ori- 
gin in Cen A as well as forecast studies of neutrino and 
gamma-ray expected signals from Cen A (e.g. Cuoco & 
Hannestad 2008, Gupta 2008, Halzen & O'Murchadha 
2008, Holder et al. 2008, Kachelriess et al. 2008). How- 
ever, such a correlation does not necessarily imply that 
Cen A is the source of these cosmic rays. First of all, 
and as noted in Abraham et al. (2008b), Hague et al. 
(2009), the evidence for clustering is based on an a pos- 
teriori analysis, so that it is difficult to quantify the level 
of significance. Even if future data confirm this cluster- 
ing, one must keep in mind that Cen A lies in front of one 
of the largest concentrations of matter in the local Uni- 
verse, at ~ 50 Mpc from us, so that the correlation could 
be accidental. Finally, as we argue in Section 14.21 the 
Cen A AGN/jet/lobe system is too weak to accelerate 
protons up to the observed energies > 57EeV. 

Even more intriguing is the fact that the chemical 
composition measured by the PAO becomes increasingly 
heavier above the ankle at 4EeV, with the last data 
points at 30 — 50 EcV lying close to expectations for iron 
(Unger et al. 2007, Wahlberg et al. 2009). Such a mea- 
surement is in obvious conflict with the most recent anal- 
ysis of HiRes data, which indicate a pure proton compo- 
sition up to 30EeV (Belz 2009). As we argue in the 
following, it is also at odds with the observed anisotropy 
pattern above 55 EeV, unless the composition abruptly 
changes from heavy to light above 30 — 50 EeV. 

3.1. Anisotropy towards Cen A 

Let us first discuss the anisotropy toward Cen A, using 
the most recent PAO dataset, in particular, 12 events lo- 
cated within 18° out of a total of 58 events above 55 EeV 
(Hague et al. 2009), assuming that the clustering is not 
a statistical accident. In a first approximation, one can 
model the anisotropy pattern as in the previous section, 
with one (or several clustered) source(s) contributing a 
fraction x of the total flux, the remaining 1 — x being ac- 
counted for by an isotropic background. We do not need 
to assume that Cen A itself is a source, but we center 
the anisotropic signal on the location of Cen A in order 
to match the observed pattern. Given that 2.7 events 
are expected on average from this region of the sky if the 
sources are distributed isotropically, we set x = 0.15 in 
the following. 

We assume that the sources inject pro- 
tons and iron nuclei with powerlaw spectra 
qz(E) oc E~ s exp{-E/[ZE majx (p)]}, producing an 
angular image of size 59. As before, this angle could 
represent smearing due to angular deflection, in which 
case 59 oc (E/Z) -1 , or the angular size of the ensemble 
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Fig. 1. — Top panel: differential histogram of the number of 
events above 55 EeV as a function of the angular separation to 
Cen A, in solid line for the model described in the text (sources 
contributing 15% of the flux, forming an image of size 88 = 10°, su- 
perimposed on an isotropic background) and expected signal from 
an isotropic population of sources (colored region). The flux is nor- 
malized to a total of 58 events above 55 EeV. Middle panel: same 
histogram but for E > 2.2 EeV, assuming that the high energy 
events are iron nuclei and that protons are injected with composi- 
tion ratio g p (E) : qz(E) = 1 : 0.06 at E < 3EeV. The anisotropy 
signal is much stronger, as measured by the number of events in 
each bin. Bottom panel: in filled circles, the ratio of the differential 
histogram shown in the middle panel to the isotropic expectations 
d7V; so /d0 with error bars estimated through the Monte Carlo; in 
empty circles, the same histogram but for a model assuming that 
the high energy events are protons only and 88 oc (E/Z)~ x . Only 
the latter model is compatible with an isotropic distribution of 
arrival directions at energies > 2.2 EeV. 



of sources. Although we assume that all heavy nuclei 
are iron, the following could be easily generalized to 
any kind of mixed composition without changing the 
basic result. We also neglect the energy losses for the 
spectrum of the anisotropic component, as before. This 
is motivated by the proximity of matter in this direction 
of the sky and, as discussed before, incorporating energy 
losses would not diminish the signal we are to extract 
out of the data, while it would introduce other free 
parameters. The remaining 1 — x fraction of the flux is 
modelled with an isotropic background, e.g. as produced 
by far away sources. 

Let us first assume that the composition is dominated 
by iron at the highest energies. We set EW = 55 EeV, 
and 59 — 10° above £W- We also set i? max (p) = 3 EeV, 
so that essentially no proton contribute to the signal 
above E t hr- In the top panel of Fig. [1] we plot the his- 
togram of the number of events as a function of the angu- 
lar distance to the source (solid line) . The colored region 
indicates the expectations for a purely isotropic signal. 
The departure from isotropy is clear, although the uncer- 
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tainty on the signal in each bin is substantial given the 
small number of events involved. Here the number of 
events within 18° of Cen A is ~ 9 for the model (slightly 
lower than the observed 12 events, yet within the uncer- 
tainty) and 2.7 for the isotropic expectations so that the 
inferred E Fc (> 55EeV) = AN/^/W~ ~ 5.5. 

In the middle panel of Fig. [TJ we plot the same his- 
togram, but at an energy E t h r /26 = 2.2EeV. In order 
to simulate the anisotropy signal, we have proceeded 
as follows. We have assumed that the total number 
of events increases (with decreasing energy) according 
to the observed data (in particular the PAO spectrum 
published in Yamamoto et al. 2007), and we have cal- 
culated the contribution to this flux above i?thr/26 for 
all the protons injected by the source above this energy. 
We have assumed that the ratio of protons to iron is 
1:0.06 at a given energy well below E max , correspond- 
ing to the composition ratio of proton to iron peak el- 
ements (Z > 17) in the galactic cosmic ray spectrum 
reconstructed at the source. The iron nuclei injected by 
the source at energies > 2.2 EeV add up to the isotropic 
background and do not contribute to the anisotropy sig- 
nal. As made obvious in Fig.[IJ the gain in signal-to- noise 
ratio is very large: the number of events observed above 
2.2 EeV within 18° is now 51 50, f or 2330 expected, so 
that S p (> 2.2 EeV) = AN/^/T^~ ~ 100. This gain by 
a factor ~ 20 can be recovered from Eq. ([2|). 

Finally, the lower panel of Fig. [1] shows the ratio of the 
differential histograms dN/d0 for the source proton con- 
tribution plus isotropic background (filled circles) to the 
isotropic background dN- mo / d9 normalized to a same to- 
tal number of events, including the error bars evaluated 
through the Monte Carlo simulation. The empty circles 
show the same quantity, but for a model in which the 
source injects only protons, even up to the highest ener- 
gies, assuming that 59 represents magnetic deflection, i.e. 
56 oc (E/Z)- 1 : we have set 50 = 10° above 55 EeV in 
this model, so that the anisotropy signal above 55 EeV is 
similar to that of the previous iron+proton model shown 
in the top panel of Fig. [IJ but the low energy signal is 
essentially isotropic, see the bottom panel. 

The conclusions to be drawn are clear. If the clustering 
seen toward Cen A is real, the events responsible for this 
anisotropy pattern toward Cen A above 55 EeV are un- 
likely to be heavy nuclei, otherwise PAO would have ob- 
served a much stronger anisotropy at energies < 40 EeV 
(Abraham et al. 2007, 2008b). 

3.2. Correlation with the Veron-Cetty & Veron (2006) 
catalog of A GN 

In the present Section, we investigate how the 
anisotropy signal with respect to nearby AGN, as re- 
ported in Hague et al. (2009) can be used to discriminate 
the chemical composition at the highest energies. We 
proceed in a similar way as before, but we now assume 
that the spectrum is composed of two components: one 
that follows the local large scale structure, contributing 
a fraction x to the flux above 55 EeV, and one isotropic 
component that makes up for the remainder, 1 — x. As a 
proxy for the local large scale structure, we use the PSCz 
catalog of galaxies (Saunders et al. 1991), as in previ- 
ous studies (Kashti & Waxman 2008, Kotera & Lemoine 
2008). Out of simplicity, we assume that the sources fol- 
low the PSCz distribution without bias. However, we 
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Fig. 2. — Differential histograms of the number of events above 
55 EeV as a function of the angular separation to the closest AGN, 
drawn from the Veron-Cetty & Veron (2006) catalog up to a dis- 
tance of 75 Mpc. The solid line corresponds to the model in which 
90% of the flux above 55 EeV is produced by sources distributed 
as the large scale structure, modeled through the PSCz catalog of 
galaxies, smeared by a deflection angle 68, as indicated. The col- 
ored region indicates the expectation for purely isotropic arrival 
directions. The histograms are normalized to a total of 58 events, 
as observed by PAO. 

allow for the possibility of a deflection angle 59 in inter- 
vening magnetic fields and we explore two possible val- 
ues 69 = 0° and 69 = 10° above 55 EeV. This choice is 
motivated by the resolution of the PSCz catalog, approx- 
imately 7°, and by the possibility that this signal is com- 
posed of heavy nuclei. In practice, we proceed through 
Monte Carlo simulations as follows. We construct many 
mock catalogs of angular sky distributions, each mock 
catalog being a sample of 58 events above 55 EeV (or 
the corresponding number at any other energy). The ar- 
rival direction of each event in the sample is drawn at 
random, according to a probability law reproducing the 
PSCz distribution and PAO coverage if it belongs to the 
x anisotropic fraction, or reproducing isotropic arrival 
directions and PAO coverage if the event belongs to the 
remaining 1 — x fraction of isotropic background. In or- 
der to draw an event from the PSCz catalog, we include 
the appropriate weight factors for the flux dilution with 
inverse distance squared, for the PSCz exposure function 
and for the PAO aperture. Then we randomly displace 
this arrival direction by 59 in the plane of the sky. The 
energy losses are accounted for by the depth of the cat- 
alog, which is limited to 200 Mpc. 

We have set x — 0.90, which allows to produce an 
anisotropy signal close to the level observed: the PAO 
indeed records 24 events out of 58 events above 55 EeV 
located within 3° of an AGN closer than 75 Mpc, while 
the above model produces 19 events within 3° if 59 = 0° 
and 14 events if 59 = 10°. The model with 90% of the 
flux above 55 EeV coming from sources distributed as the 
large scale structure with negligible deflection thus gives 
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Fig. 3. — Ratio of the differential histograms for the models 
shown in Fig.f3]to histograms for purely isotropic arrival directions, 
but computed from the proton component at energies > 2.2 EeV. 
The uncertainties measured through the Monte Carlo are indicated. 
Note the change of scale between the two panels corresponding to 
the two values of 88, as indicated. 



a number of correlating AGN within 3° that is compa- 
rable to that seen by PAO, to within the uncertainty. 
A deflection angle 86 = 10° provides a signal that is 
marginally too low, but given the simplifying assump- 
tions made above, this is not crucial. It certainly indi- 
cates that larger deflection angles will probably result in 
too small anisotropics. It is interesting to note that the 
anisotropy signal peaks at a separation which is substan- 
tially smaller than 56, i.e. ~ 4 - 5° for 86 = 10°. This 
of course results from the PAO procedure, which selects 
the closest AGN to the arrival direction and not to the 
source direction, i.e. after the direction at the source has 
been modified by 86. The differential histograms corre- 
sponding to these models and isotropic expectations are 
shown in Fig. [2] as a function of the angular separation to 
the closest AGN drawn from the Veron-Cetty & Veron 
(2006) catalog, with an arbitrary cut at 75 Mpc, follow- 
ing the PAO analysis. These histograms are shown for 
the two possible values of 86, as indicated. Here, we find 
£fo(> 55 EeV) ~ 2.1 and 0.76 for 86 = 0° and 86 = 10° 
respectively. The signal is thus marginal in terms of de- 
parture from isotropy. 

We now plot the ratios of the differential histograms 
expected in these two models at low energy > 2.2 eV 
to the histogram expected for purely isotropic arrival 
directions, as before. We have assumed a conserva- 
tive composition ratio 1:0.06 between proton and iron 
as in the previous section, and we have neglected energy 
losses on the propagated spectra. The results are shown 
in Fig. [31 which clearly reveals a stronger anisotropy 
signal at these energies. In terms of signal-to-noise, 
we find S p (> 2.2 EeV) ~ 6.2 and 2.3 respectively for 
86 = 0° and 86 = 10° (AN = 623 and 235 respectively, 
Ajso = 10100), which significantly exceed the signal-to- 



noise ratios obtained above 55 EeV. This generalizes the 
argument presented in Section [2] for the case of a more 
complex distribution of point sources on the sky, su- 
perimposed on an isotropic background. It also shows 
that the present angular data of the PAO disfavors a 
heavy composition at ultra-high energies. Strictly speak- 
ing, this test demonstrates that the anisotropic part of 
the signal cannot be due to heavy nuclei but one can- 
not exclude that the isotropic background is made up of 
heavy nuclei. However, our above model also suggests 
that a significant fraction (90% in Fig. I2I3| of the flux 
at > 55 EeV is produced by the anisotropic component if 
one takes at face value the anisotropy results of the PAO. 
This suggests that the bulk of ultra-high energy cosmic 
rays above 55 EeV are protons. 

4. AGN AS SOURCES OF ULTRA-HIGH ENERGY COSMIC 

RAYS? 

4.1. AGN luminosity vs acceleration 

On general grounds, one can construct a lower bound 
to the magnetic luminosity that a source must exhibit 
in order to be able to accelerate particles up to 10 20 eV 
(Norman et al. 1995, Waxman 1995, Lyutikov & Ouyed 
2005, Waxman 2005). Let us first briefly describe the 
general, model independent argument given in Waxman 
(1995). Consider an astrophysical source driving a flow 
of magnetized plasma, with characteristic magnetic field 
strength B and velocity v. Imagine now a conducting 
wire encircling the source at radius R. The potential 
generated by the moving plasma is given by the time 
derivative of the magnetic flux, and is therefore given by 
V = /3BR where (3 — v/c. A proton which is allowed 
to be accelerated by this potential drop would reach 
energy E p ~ (3eBR. The situation is somewhat more 
complicated in the case of a relativistic outflow, where 
r = (1 — ft 2 )^ 1 / 2 ^> 1. In this case, the proton is allowed 
to be accelerated only over a fraction of the radius R, 
comparable to R/T. To see this, one must realize that as 
the plasma expands, its magnetic field decreases, so the 
time available for acceleration corresponds, say, to the 
time of expansion from R to 2R. In the observer frame 
this time is R/c, while in the plasma rest frame it is R/Tc. 
Thus, a proton moving with the magnetized plasma can 
be accelerated over a transverse distance ~ R/T. This 
sets a lower limit to the product of the magnetic field 
and source size, which is required to allow acceleration 
to E p , BR > TEp/e/3. This also sets a lower limit to 
the rate Lb at which energy is carried by the out flow- 
ing plasma, and which must be provided by the source 
(Waxman 1995, 2005), 

r 2 / 171 \ 2 r2 / ,ry \ 2 



L B > 



E 
~Ze 



10 



45.5 1 



E/Z 

10 20 eV 



erg/s, (3) 



P \ZeJ 

where we have generalized the last equation to particles 
of charge Z. 

Let us consider next the application of the general ar- 
gument given above to the commonly assumed scenario, 
where particle acceleration within the outflowing plasma 
is achieved via diffusive acceleration. In this case, one 
writes the acceleration timescale t acc = At^ as a mul- 
tiple A of the Larmor time, and one compares it to the 
dynamical timescale available for acceleration. Let us 
discuss this bound in some detail. For the sake of defi- 
niteness, we assume an outflow of Lorentz factor T and 
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opening angle (the isotropic case is recovered in the 
limit 2 — ► 2). The quantity A depends on the diffu- 
sion coefficient of the particle, hence on its energy and 
the general properties of the magnetic field of the ac- 
celerating region. In acceleration mechanisms which in- 
volve scattering against magnetic inhomogeneities (such 
as Fermi I or Fermi II), i acc oc t s , where t s = 3k/ c 
is the scattering timescale and k the diffusion coeffi- 
cient. The parallel diffusion coefficient can be written 
as k/(t"lc) ~ r/ -1 (r^/As)" (see for instance Casse, 
Lemoine & Pelletier 2002), with -q < 1 the fraction of 
turbulent magnetic energy density (in units of total mag- 
netic energy density) and Xb the coherence length. The 
exponent a depends on the spectral index (3 of the (one 
dimensional) turbulence power spectrum: a = (3 — 1 for 
t-l < Xb, a = 1 for tl > Xb- From this, one easily 
verifies that k/(t"lc) > 1 is general, meaning that the 
so-called Bohm regime k = r-^c is a limiting regime. In 
non-relativistic Fermi shock acceleration (r s h/3 s h < 1 or 
Ah ^5 0.7), t acc ~ 3«;//3g h , with /3 s h the shock velocity in 
the comoving frame and k refers to quantities measured 
in the upstream frame. Hence A ~ 3/3~^ k/ '(jic) ^> 1 
. One recovers a similar result for non-relativistic Fermi 
II acceleration, with A <~ f3~^ k/{tlc) 3> 1, Pa < 1 
denoting the Alfven speed. If acceleration takes place 
in a (non-relativistic) shear flow with velocity gradi- 
ent Am/ Ax, the timescale i acc <~ Ax 2 j (Au 2 t s ) (Rieger, 
Bosch- Ramon & Duffy 2007). The deconfinemcnt limit 
corresponds to t a ~ Ai/c, in which case the limiting 
acceleration timescale becomes comparable to t s c 2 /Au 2 , 
as above. Note that the acceleration timescale is larger 
(hence acceleration is less efficient) when t s < Ax/c. In 
the ultra-relativistic shock limit r s h/3 s h S> 10, one ex- 
pects A ~ T^T^/As (up to a factor of order unity, see 
Pelletier, Lemoine & Marcowith 2009), with the condi- 
tion tl/A_b 3> r s h- Since A increases with n,, this pro- 
cess is in general not efficient. Finally, in the mildly rela- 
tivistic regime 1 < r s h/3 s h ^ 10, one expects t acc ~ t s , 
hence A ~ k/{tlc), whose minimum value is of order 
one provided the Bohm regime diffusion applies and the 
magnetic field is nearly fully turbulent. 

To summarize, A 1 in general, and A « 1 cor- 
responds to a maximally efficient acceleration process. 
The constraint i acc < tdyn, with tdyn = R/ (3Tc the dy- 
namical timescale at radius R (in the lab frame), now 
gives: 

L B > 0.65 x 10 45 Q 2 Y 2 A 2 f Z~ 2 E 2 2Q erg/s , (4) 

with E20 the observed energy in units of 10 20 eV. The 
presence of /3 3 does not mean that in the limit (3 — > 0, 
this lower bound vanishes, since A 2 oc [3~^ more than 
compensates for this term. Similarly, the above bound 
does not vanish in the limit — > 0, since lateral es- 
cape losses become prominent when O < The 
corresponding timescale t csc ~ (OR) 2 /2k, hence Lb > 
1.2 x 10 45 ^/3(K/r L c)Z- 2 £;| crg/s. In the case of AGN 
jets, one finds in general Or ~ 1, hence both limits are 
comparable. According to the above discussion, the most 
optimistic values for the parameters entering these equa- 
tions are /3 ~ 1, A ~ 1 and k/(t-lc) <~ 1. Therefore 
10 45 Z~ 2 erg/s can be considered as a firm lower bound 
on the source magnetic luminosity. 

Since most of the AGN seen to correlate within 3° of 
the PAO events are Seyfert galaxies, with bolometric lu- 



minosities well below 10 ergs/s, the AGN that they har- 
bor cannot be the source of light particles at 10 20 eV. We 
also note that Zaw, Farrar & Greene (2009) have shown 
that one third of this sample is actually made of star 
forming galaxies with little or no AGN activity. Hence 
the correlation with stricto-sensu AGN objects is actu- 
ally weaker than reported. 

Actually even FR I radio-galaxies, TeV blazars and BL 
Lac objects do not possess the required power to accel- 
erate Z ~ 1 particles up to 10 20 eV, since their magnetic 
luminosities are of order Lb ~ 10 42 — 10 44 ergs/s (Celotti 
& Ghiscllini 2008). Such objects might possibly acceler- 
ate heavy nuclei up to energies around the GZK cut-off 
but, as discussed above, the anisotropy pattern at high 
and especially low energies would then appear inconsis- 
tent with the PAO results. Further anisotropy studies 
in both energy ranges would certainly further constrain 
their possible contribution. 

The compilation of Celotti & Ghisellini (2008) has been 
done in the frame of leptonic synchrotron and inverse 
Compton emission. In hadronic blazar models, the mag- 
netic field in the blazar zone is typically one order of 
magnitude larger than in leptonic models, so that these 
objects might fulfill the magnetic luminosity bound to 
accelerate protons up to 10 20 eV. However, in this case, 
acceleration would occur in the blazar zone hence the 
emission should be beamed forward. In order to escape 
further expansion losses in the magnetized jets, the ac- 
celerated protons would have to be converted into neu- 
trons, which would decay back to protons on a distance 
scale ~ O.9E20 Mpc, i.e. outside the jet. One should 
therefore observe a correlation of the arrival directions 
with blazars, not with radio-galaxies seen offside (Rachcn 
2008). However the Pierre Auger Observatory reports no 
correlation with blazars (Harari et al. 2007), so that this 
possibility is equally at odds with present observational 
results. 

Keeping in line with the modelling of Celotti & Ghis- 
cllini (2008), only flat spectrum radio quasars (the likely 
equivalent to the most powerful FR II sources) seem ca- 
pable of producing jets with Lb > 10 45 ergs/s. But, 
considering the highest energy PAO event (E = 1.48 ± 
0.27 ± 0.32 x 10 20 eV), one finds that the smallest angu- 
lar separation of this event to the FR II sources located 
within 130 Mpc, as compiled by Massaglia (2007) is al- 
ready 28° (for NGC 4261). The next objects are 3C 296 
(38° away) and PKS1343-60 (41° away). The closest 
blazar (classified as BL Lac in the Veron & Veron-Cetty 
catalog) located closer than 150 Mpc lies 115° away from 
this highest energy event (TEX 0554+534). 

Recently, it has been proposed that transients in active 
galactic nuclei could power up the engine up to the lumi- 
nosities required, thereby evading the above constraints 
(see for instance Farrar & Gruzinov 2009, Dermer et al. 
2008). However, as demonstrated by Waxman & Loeb 
(2009), such transients should produce counterpart flares 
in X-rays through the concomittant acceleration of elec- 
trons, which should have been picked up by existing sur- 
veys. Their non-detection strongly argues against such 
flaring scenarios. 

4.2. The particular case of Cen A 

Given its proximity and the apparent clustering of 
PAO events in this direction of the sky, Cen A has re- 
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ceived a lot of attention. In particular, several authors 
have suggested that this AGN could host a site of acceler- 
ation to ultra-high energies (Cavallo 1978; Romero et al. 
1996; Farrar & Piran 2000; Gorbunov et al. 2008; Hard- 
castle et al. 2009, O'Sullivan, Reville & Taylor 2009). 
Let us therefore discuss it in some more detail. 

Centaurus A is classified as a FR I / misaligned blazar 
(see Israel 1998 for a review). Its bolometric luminos- 
ity L bo \ ~ 10 43 erg/s and its jet kinetic power is in- 
ferred to be Lj e t ~ 2 x 10 43 erg/s. Therefore, this source 
should not be able to accelerate light particles to 10 20 eV. 
The detailed modelling of the spectral energy distri- 
bution of its nucleus by Chiaberge, Capetti & Celotti 
(2001) provides further information on the emission zone: 
R ~ 10 16 cm, B ~ 0.5 G and V < 3 - 5. This corre- 
sponds to a magnetic luminosity Lb = \R 2 B 2 T 2 j3c ~ 
2 x 10 42 ergs/s, a value that seems quite reasonable in 
view of the jet kinetic power and of the bolometric lu- 
minosity, and which agrees relatively well with estimates 
of the magnetic field in the inner jet at larger distances 
from the core, when assuming BR ~ constant. 

The papers which have discussed the issue of accelera- 
tion to ultra-high energy in some detail are Romero et al. 
(1996), which has argued that acceleration can take place 
in the X-ray knots of the inner jet through diffuse shock 
acceleration, Hardcastle et al. (2009) and O'Sullivan, Re- 
ville & Taylor (2009), which have discussed the stochastic 
acceleration of protons in the giant lobes of Ccn A and 
Gorbunov et al. (2008) relying on the model of Neronov 
et al. (2007) of particle acceleration in a black hole mag- 
netosphere. 

Romero et al. (1996) uses the inferred values of 
B ~ 10~ 5 G (upstream field), R ~ 1.8 kpc, /3 sh ^ 
0.3 and r\ ~ 0.4. However, since the Larmor radius 
r L ~ 11\ivcE 2 qZ- 1 BzI > R (with E 2Q = £/10 20 eV, 
S_ 5 = B/10- 5 G), even the Hillas criterion (Hillas 1984) 
is not satisfied. Balancing the acceleration timescale 
with the escape timescale R/c, one derives a maximal 
energy E < 10 18 Z0 1 ? B_ 5 i?i. 8 cV {R x .% = i?/ 1.8 kpc, 
/3o.3 = /3 s h/0.3). As discussed above, this maximal en- 
ergy lies in the correct ballpark for accelerating nuclei to 
ultra- high energies. However the anisotropy signal would 
then appear incompatible with the PAO data; in partic- 
ular, a strong anisotropy signal at EeV energies should 
have been detected. 

Hardcastle et al. (2009) have argued that protons 
could be accelerated up to <~ 5 x 10 19 eV through stochas- 
tic acceleration in the magnetized turbulence of the gi- 
ant lobes of Cen A. However, one can show that their 
estimate of the acceleration timescale is far too opti- 
mistic because, when applied to the electrons, it pro- 
duces a maximal energy well in excess of the inferred 
£ maX]e <~ 1 — 4 x 10 11 eV. In detail, using the esti- 
mate of Hardcastle et al. (2009) for the acceleration 

timescale, i acc ~ 3.5 x 10 6 yrs^o^Zgj and their Bohm 
scaling i acc cx E, one infers the maximal energy for elec- 
trons by comparing i acc to the synchrotron and inverse 
Compton energy loss timescales, leading to -E max ,e ~ 

2 x 10 16 B_y 2 cV. The above maximal electron energy is 
five orders of magnitude larger than the inferred maxi- 
mal energy, hence the above acceleration timescale is far 
too optimistic and, consequently, the maximal proton en- 
ergy must be much lower than 10 20 eV. This conclusion 



agrees with those of O'Sullivan, Reville & Taylor (2009); 
see also Casse, Lemoine & Pelletier (2002) for a simi- 
lar discussion regarding the acceleration in the lobes of 
powerful radio-galaxies. 

Finally, Neronov et al. (2007) have argued that in a 
magnetic field of strength B = 10 4 G threading a maxi- 
mally rotating black hole of mass Mbh = 10 8 M Q , heavy 
nuclei could be accelerated to energies of order 10 20 eV 
by exploiting the potential drop $ ~ 10 20 i?4M 8 V. How- 
ever, this maximal energy is quite optimistic because it 
ignores radiative energy losses, which are highly efficient 
in the nuclei of powerful AGN (Norman et al. 1995, 
Henri et al. 1999). In any case, the Poynting flux that 
is required, Lb ~ 10 45 erg/s far exceeds any other lu- 
minosity measured elsewhere in the galaxy, and in par- 
ticular the estimate obtained from multi-band spectral 
modelling of the nucleus emission (Chiaberge, Capetti & 
Celotti 2001). 

5. DISCUSSION 
5.1. Summary 

Let us start by summarizing the results obtained so 
far. We have discussed how one can can test the chemi- 
cal composition of ultra-high rays on the sky by compar- 
ing the anisotropy signals at various energies. Our main 
result is to show that if anisotropies are observed above 
some energy -Ethr and the composition is assumed to be 
heavy at that energy (nuclei of charge Z), one should 
observe at energies > E^/Z a substantially stronger 
anisotropy signal, which is associated with the proton 
component emitted by the sources that are responsible 
for the anisotropy pattern. These conclusions do not de- 
pend at all on the modelling of the intervening magnetic 
fields and they are robust with respect to the parameters 
characterizing the sources. 

We have then applied this test to the most recent data 
of the PAO. We find that, if these data are taken at face 
value (and if the PAO arrival directions around the ankle 
appear isotropic) , the events that are responsible for the 
anisotropy signal reported toward Cen A and toward the 
nearby AGN should not be heavy nuclei. Interestingly, 
this result appears to be at odds with the current PAO 
results on the chemical composition at the highest ener- 
gies. One cannot exclude at the present time that the 
observed anisotropies are statistical accidents; only the 
acquisition of a larger set of data will tell. One cannot 
exclude either that the composition switches abruptly 
from heavy to light at <~ 30EcV. Finally, one cannot 
exclude that a systematic bias affects the composition 
measurements of the PAO. In any case, the present dis- 
cussion indicates that, if both the anisotropy and the 
heavy composition are confirmed by future data, some 
possibly important information is to be extracted out of 
this apparent contradiction. 

The theoretical and experimental uncertainties in the 
extrapolation of the p — p cross-section to center-of-mass 
energies y/s > 100 TeV are one of the possible sources 
of biases in shower reconstruction. Recently, it has been 
argued that, if this cross-section were underestimated 
by some 40 — 60 % at these energies, one might recon- 
cile the existing A max measurements with a pure proton 
composition at energies above the ankle (Ulrich et al. 
2009a, 2009b; see however Wibig 2009). In this respect, 
and in the light of the above discussion, the apparent 
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anisotropy towards Cen A offers a very valuable oppor- 
tunity to test this possibility, ft suffices, in principle, 
to test the anisotropy pattern against a future (indepen- 
dent) data set at various energies by comparing the num- 
ber of events detected within a predetermined solid angle 
in a predetermined direction to the expected number for 
isotropic arrival directions; at the same time, one needs 
to perform a dedicated composition measurement for cos- 
mic rays from this particular region of the sky, in order 
to avoid as much as possible the contamination due to 
the isotropic contribution to the all-sky flux. Then, one 
should apply the test discussed above. 

5.2. The interpretation of the PAO data 

The natural question that follows is what can be in- 
ferred about the sources of ultra-high energy cosmic rays 
if, following the anisotropy results of PAO and the above 
discussion, one assumes that most ultra-high energy cos- 
mic rays are protons. Let us now explore this line of 
reasoning. 

First of all, one concludes that the existing data of the 
PAO disfavor the AGN model of ultra-high energy cos- 
mic ray origin. Indeed, we have discussed in Section |4] 
a robust lower bound to the magnetic luminosity that a 
source must possess in order to be able to accelerate par- 
ticles of charge Z up to 10 20 eV, L B > 1 45 Z~ 2 erg/s. 
This bound allows to rule out the local radio quiet AGN 
as sources of 10 20 eV protons. We have also argued that, 
while FR I radio-galaxies do not seem to possess the re- 
quired luminosity to accelerate protons up to I0 20 eV, 
FR II radio-galaxies and blazars do not correlate with 
the highest energy events. Although radio-loud galax- 
ies could possibly accelerate heavy nuclei up to energies 
close to the GZK cut-off, as we have discussed they would 
produce an anisotropy pattern predominantly at ~ EeV 
energies rather than at ultra-high energies, in conflict 
with the PAO results. 

We have discussed the particular case of Cen A in some 
detail, as this source has recently received a lot of at- 
tention due to the observed clustering of events around 
it. We have provided a critical discussion of acceleration 
in this object and argued that its AGN/jet/lobe system 
cannot produce protons at the energies required, in good 
agreement with the general expectations from the above 
luminosity bound. That being said, the excess clustering 
in this particular direction, if real, remains to be inter- 
preted. 

First, if sources are distributed according to the large 
scale structure, one expects a certain number of events 
to coincide accidentally with the direction of Cen A, be- 
cause of the location of the Centaurus and Shapley su- 
perclusters in this direction. Using the PSCz catalog of 
galaxies as a proxy for the distribution of the sources, 
with a depth of 200 Mpc, and accounting for the PAO 
exposure, one expects 0.80 events above 55 EeV within 
6° of Cen A nucleus (and 0.3 for isotropic arrival direc- 
tions) for 2 — 3 observed, or 4.4 events within 18° (and 
2.7 for isotropic arrival directions) for 12 observed. Note 
that the large scale structure is poorly sampled in the 
proximity of Cen A due to its low galactic latitude. Nev- 
ertheless, the above numbers already indicate that the 
contamination from the surrounding large scale structure 
in the direction to Cen A is substantial. In particular, 
at the present level of statistics, one cannot claim ob- 



serving a significant departure from a model in which 
the sources are distributed as the large scale structure. 
It should also be recalled here that the PAO analysis of 
clustering around Cen A is an a posteriori one, hence 
one may not therefore assign a reliable significance to 
the detection of the anisotropy. 

Moreover, if sources are gamma-ray bursts or magne- 
tars, one should include in the flux prediction the contri- 
bution of such sources located in the Cen A host galaxy 
itself. It is well known that the probability of detect- 
ing ultra-high energy cosmic rays from nearby gamma- 
ray bursts is extremely small unless the arrival times 
of these ultra-high energy cosmic rays are sufficiently 
dispersed (Waxman 1995, Waxman & Miralda-Escude 
1996). Consider for instance gamma-ray bursts with rate 
of occurence Nqbb in Cen A; these can be seen if ei- 
ther of their jets points into a solid angle 2ir(l — cos 89) 
where 89 is the typical deflection acquired by the cosmic 
ray en route to the detector, giving a detection probabil- 
ityP ~ N GBB a t S9 2 /2 ~ IO-^ghb.-s^to S^A?^, 
with N GRB = 10 _5 AGRB,-5yr~ 1 ; cr t = 6 x 
10 3 yri3 2 8 Aiookpc-£'7o 2 represents the spread of arrival 
times and 89 = Z° B -^X^^. E 70 the angular de- 
flection for particles of energy E = 70 Ero EeV trav- 
eling from Cen A to the detector through a mag- 
netic field of strength 10 _8 S_8 G and coherence length 
100 Aiookpckpc (Waxman & Miralda-Escude 1996). The 
above probability suggests that one should not detect 
events from gamma-ray bursts emitted directly toward 
the observed. However, if the jet of the gamma-ray burst 
hits one of the lobes of Cen A, the particles can be redis- 
tributed in all directions, in particular towards the de- 
tector. The magnetic field inside the lobe, B ~ 1 /iG is 
indeed sufficient to impart a deflection of order unity on 
the lobe distance scale ~ 100 kpc to protons of 70 EeV. 
In this case, the expected time delay and time disper- 
sion are of the order of the travel time to the lobes, 
L\ohes/c ~ 3 x 10 5 yr, so that the expected number of 
gamma-ray bursts that one can see at any time, through 
rescattering of the particles on the lobes of Cen A is: 

(V GRB ) ~ V GRB Ll ° bcsA " lobcs ~ 0(1)Vgrb,- 5 (5) 

C 27T 

for an apparent solid angle of the lobes Af2i o b 0S as viewed 
from the host galaxy of Cen A, which is close to Air. 

The last gamma-ray burst (s) of Cen A may then con- 
tribute to the flux above the GZK cut-off up to: 

JCenA « 0.9x KT 41 eV-WV 1 esi^ro/fVe 5 (N GRB ) , 

(6) 

with er ti5 .5 = Of/3 x 10 5 yrs, f e = ln(E max /E m i n ) ~ 
1 — 10 for an injection spectrum with s — 2.0, and 
e = 10 51 £5i ergs the energy injected by one gamma-ray 
burst in ultra-high energy cosmic rays. This flux should 
be spread over the angular scale of the lobes, ~ 10°. 
At 70 EeV, PAO records a diffuse flux j PAO (70EeV) ~ 
3 x 10~ 40 eV cm - 2 8 _1 sr -1 . After correcting for the 
solid angle of the Cen A image and the PAO aperture to- 
wards Cen A, one finds that the last gamma-ray bursts 
in Cen A should contribute up to 2 — 25 % of the ob- 
served PAO flux, well within the range of the observed 
anisotropy. 

We also note that this contribution from sources in- 
side Cen A would improve the overall correlation of all 
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the PAO data above 55 EeV with sources distributed as 
the local large scale structure. Finally, we remark that 
the above effect of rescattering on the lobes provides a 
clear example of the possible biases discussed in Kotera & 
Lemoine (2008), namely that the PAO could be preferen- 



Abbasi, R. U., et al. (The HiRes Collaboration), 2004, ApJ, 610, 
L73 

Abbasi, R. U., et al. (The HiRes Collaboration), 2005, ApJ, 623, 
164 

Abbasi, R. U., et al. (The HiRes Collaboration), 2008, Physical 

Review Letters, 100, 101101 
Abraham, J. et al. (The Pierre Auger Collaboration), 2007, Science, 

318, 938 

Abraham, J., et al. (The Pierre Auger Collaboration), 2008a, 

Physical Review Letters, 101, 061101 
Abraham, J. et al. (The Pierre Auger Collaboration), 2008b, 

Astroparticle Physics, 29, 188 
Allard, D., Busca, N. C, Decerprit, C, Olinto, A. V., & Parizot, 

E. 2008, Journal of Cosmology and Astro-Particle Physics, 10, 

33 

Allard, D., & Protheroe, R. J. 2009. larXiv:0902.4538l 

Aloisio, R., Berezinsky, V., & Gazizov, A. 2008, arXiv: 0803. 2494 
Arons, J. 2003, ApJ, 589, 871 

Belz, J. (The HiRes Collaboration), proc. CRIS-2008 (Malfa, Italy, 
Sept. 2008) 

Aublin , J. et a l. (The Pierre Auger Collaboration), 2009, 

larXiv:0906.2347 ICRC 2009 
Bertone, G., Isola, C, Lemoine, M., & Sigl, G. 2002, Phys. Rev. D, 

66, 103003 

Casse, F., Lemoine, M., & Pelletier, G. 2002, Phys. Rev. D, 65, 
023002 

Cavallo, G. 1978, A&A, 65, 415 

Celotti, A., Ghisellini, G., 2008, MNRAS, 385, 283 

Chiaberge, M., Capetti, A., & Celotti, A. 2001, MNRAS, 324, L33 

Cuoco, A., & Hannestad, S. 2008, Phys. Rev. D, 78, 023007 

D ermer, C. P., Ra zzaque, S., Finke, J. D., & Atoyan, A., 2008, 

larXiv:0811. 11601 
Evans, N. W., Ferrer, F., Sarkar, S., 2003, PRD 69, 3501. 
Evans, N. W., Ferrer, F., Sarkar, S., 2004 , PRD 69, 8302. 
Farrar, G. R., & Piran, T. 2000, |arXiv:astro-ph/0010370| 
Farrar, G. R., Berlind, A. A., & Hogg, D. W. 2006, ApJ, 642, L89 
Farrar, G. R., & Gruzinov, A., 2009, ApJ, 693, 329 
Finley, C. B., & Westerhoff, S. 2004, Astroparticle Physics, 21, 359 
Ghisellini, G., Ghirlanda, G., Tavecchio, F., Fraternali, F., & 

Pareschi, G. 2008, arXiv:0806.2393 
Gorbunov, D. S., Tinyakov, P. G., Tkachev, I. I., Troitsky, S. V., 

2002, ApJ 577, L93. 
Gorbun ov, D. S., Tinyakov , P. G., Tkachev, I. I., Troitsky, S. V., 

2004, astro-ph/0406654 
Gorbun ovTD. S., Tinyako v, P. G., Tkachev, I. I., & Troitsky, S. V. 

2008, arXiv:0804.1088 
Greisen, K., 1966, PRL 16, 748. 

Gupta, N. 2008, Journal of Cosmology and Astro-Particle Physics, 
6, 22 

Hague, J. D. et al. (The Pierre Auger Collaboration), 

larXiv:0906.2347 ICRC 2009 

Halzen, F., & O'Murchadha, A. 2008, arXiv:0802.0887 

H arari, P., et al . (The Pierre Auger Collaboration), 2007, 

larXiv:0706.1715l 
Hardcastle, M. J., Cheung, C. C, Feain, I. J., & Stawarz, L. 2009, 

MNRAS, 393, 1041 
Henri, G., Pelletier, G., Petrucci, P. O., & Renaud, N. 1999, 

Astroparticle Physics, 11, 347 
Hillas, A. M. 1984, ARA&A, 22, 425 

H older, J., fe for the VERITAS Collaboration 2008, 
larXiv:0810.047T1 



tially observing the last scattering centers (in the present 
case, the lobes of Cen A) rather than the source itself (the 
gamma-ray bursts located in the core of Cen A). 



Israel, F. P. 1998, A&A Rev., 8, 237 

Kachelriess, M., Ostapchenko, S., & Tomas, R. 2008, 

[^Civ:0805.260|] 
Kashti, T., & Waxman, E. 2008, Journal of Cosmology and Astro- 
Particle Physics, 5, 6 
Kotera, K., & Lemoine, M. 2008, Phys. Rev. D, 77, 1 23003 
Lyutikov, M., & Ouyed, R. 2005, |arXiv7ast ro-ph/0507620 
Massaglia, S., 2007, Nuc. Phys. Bl Proc. Supp. 168, W2 
Milgrom, M., & Usov, V. 1995, ApJ, 449, L37 

Moskal enko, I. V., Stawa rz, L., Porter, T. A., & Cheung, C. C. 

2008. [arXiv:0805. 12601 

Nagar, N. M., & Matulich, J. 2008. larXiv:0806.3220l 

Neronov, A., Semikoz, P., & Tkachev, I. 2007, |arXiv:0712.1737 
Norman, C. A., Melrose, D. B., & Achterberg, A. 1995, ApJ, 454, 

60 

O'Sullivan, S., Reville, B., & Taylor, A. M., 2009, arXiv:0903.1259 
Pelletier, G., Lemoine, M., & Marcowith, A., 2008, arXiv:0707.3459 
P anov, A. P., et al. (ATIC-2 Collaboration) 2006, 

arXiv:astro-ph/0612377 
Rachen, J. P., & Biermann, P. L. 1993, A&A, 272, 161 
Rachen, J., 2008, in Proceedings of "Rencontres dc Blois 2008 - 

Challenges in particle astrophysics", Blois (France), May 18-23 

2008, ed. J. Pumarchez 
Rieger, F. M., Bosch-Ramon, V., & Puffy, P. 2007, Ap&SS, 309, 

119 

Romero, G. E., Combi, J. A., Perez Bergliaffa, S. E., & 
Anchordoqui, L. A. 1996, Astroparticle Physics, 5, 279 

Saunders, W., Frenk, C, Rowan-Robinson, M., Lawrence, A., & 
Efstathiou, G. 1991, Nature, 349, 32 

Stanev, T., Biermann, P. L., LLoyd- Evans, J., Rachen, J. P., 
Watson, A. A., 1995, PRL 75, 3056. 

Takeda, M., et al. 1999, ApJ, 522, 225 

Tinyakov, P. G., Tkachev, I. I., 2001, Pisma Zh. Eksp. Teor. Fiz. 

74, 3 [JETP Lett. 74, 1]. 
Tinyakov, P. G., Tkachev, I. I., 2002, Astropart. Phys. 18, 165. 
Tinyakov, P. G., Tkachev, I. I., 2004, PRP 69, 8301. 
Uchihori, Y., Nagano, M., Takeda, M., Teshima, M., LLoyd-Evans, 

J., Watson, A. A., 2000, Astropart. Phys. 13, 151. 
Ulrich, R., Engel, R., Miiller, S., Pierog, T., Schiissler, F., & Unger, 

M. 2009a, arXiv:0906.0418 
U lrich, R., Engel, R ., Miiller, S., Schiissler, F., & Unger, M. 2009b, 

[arXiv:0906. 30751 
U nger, M. et al. ( The Pierre Auger Collaboration), ICRC-2007, 

[arXiv:0706.1495] 
Veron-Cetty, M.-P., & Veron, P. 2006, A&A, 455, 773 
Vietri, M. 1995, ApJ, 453, 883 

Wahlberg, H. et al. (The Pierre Auger Collaboration), 2009, 

larXiv:0906.2319 ICRC 2009 
Waxman, E. 1995, Phys. Rev. Lett., 75, 386. 
Waxman, E., & Miralda-Escude, J. 1996, ApJ, 472, L89 
Waxman, E., 2005, Phys. Script a, T121, 147 
Waxman, E., & Loeb, A., 2008, arXiv:0809.3788 
Wibig, T. 2008, arXiv:0810.5281 
Wibig, T. 2009, Phys. Rev. P, 79, 094008 

Yamamoto, T. et al. (The Pierre Auger Collaboration), ICRC 2007 
Zatsepin, G. T., Kuzmin, V. A., 1966, Pis'ma Zh. Eksp. Teor. Fiz. 

4, 114 [JETP. Lett. 4, 78]. 
Zaw, I., Farrar, G. R., & Greene, J. E. 2009, ApJ, 696, 1218 



